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Abstract. We investigate mixing effects on the glass state of binary colloidal hard-sphere-like mixtures with large size 
asymmetry, at a constant volume fraction = 0.61. The structure, dynamics and viscoelastic response as a function of mixing 
ratio reflect a transition between caging by one or the other component. The strongest effect of mixing is observed in systems 
dominated by caging of the large component. The possibility to pack a large number of small spheres in the free volume left by 
the large ones induces a pronounced deformation of the cage of the large spheres, which become increasingly delocalised. This 
results in faster dynamics and a strong reduction of the elastic modulus. When the relative volume fraction of small spheres 
exceeds that of large spheres, the small particles start to form their own cages, slowing down the dynamics and increasing the 
elastic modulus of the system. The large spheres become the minority and act as an impurity in the ordering beyond the first 
neighbour shell, i.e. the cage, and do not directly affect the particle organisation on the cage level. In such a system, when 
shear at constant rate is applied, melting of the glass is observed due to facilitated out-of-cage diffusion which is associated 
with structural anisotropy induced by shear. 
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INTRODUCTION 

Many different systems, among them polymers, metals 
and colloids, can form thermodynamically equilibrated 
states, but also non-equilibrium, metastable states, in- 
cluding amorphous solid materials called glasses [1, 2]. 
The glass transition is generally associated with a dra- 
matic slowing down of the particle dynamics which is 
driven by changes in thermal energy or crowding. 
One of the simplest model systems to study crowding in- 
duced glass formation are suspensions of colloidal hard 
spheres. By increasing the particle volume fraction <p for- 
mation of a glass state above (j) — (j) g prevents crystallisa- 
tion, if the system has a sufficiently broad distribution 
of sizes. The formation of the glass state is explained 
in terms of the cage effect: At <p > <p g each particle is 
trapped in the cage of its neighbours resulting in dynam- 
ical arrest, i.e. the absence of long distance diffusion over 
a large window of times [3, 4]. Dynamical arrest and for- 
mation of a solid state above > ty g are also manifested 
in the viscoelastic properties as a sudden increase of the 
viscosity [5] and the appearance of a Maxwell plateau 
modulus in the linear response [6]. 

The addition of a second component with a signifi- 
cantly different mean size compared to the first compo- 
nent, leads to an even richer scenario. Depending on the 
total volume fraction of the system and the mixing ratio 
of the two species, mode coupling theory (MCT) predicts 



the existence of different glass states [7]. When the size- 
ratio 8 = R s /Ri = d s /di, where R s , d s and R/, d[ are the 
radii and diameters of the small and large components re- 
spectively, becomes about 0.2 and smaller, four different 
glass states are expected [7]: In the first state both com- 
ponents are caged; in the second state dynamical arrest 
of the large component is driven by depletion attraction 
induced by the small species; in the third state the large 
component is arrested through caging, while the small 
component is mobile; finally the small particles can be 
caged, while the large particles are not caged, but only 
localised by the surrounding dense matrix of small parti- 
cles. Despite the rich behaviour predicted by theory, the 
glass state of colloidal binary mixtures at such large size 
disparities is hardly studied experimentally [8]. In [8] the 
formation of a glass despite the mobility of the small 
component is reported. A similar glass state has also been 
found in simulations of soft sphere mixtures [9, 10]. 

In order to extend these studies and to explore the for- 
mation of different glasses, we performed experiments to 
determine the microscopic structure, dynamics and vis- 
coelastic response of colloidal hard-sphere mixtures of 
large size disparity (5 = 0.2) and constant total volume 
fraction ~ 0.61. We vary the relative volume fraction 
of the small component, x s = </>. v /0, to explore the effect 
of mixing on the glass state. We find that the composition 
of the mixture strongly affects the dynamics and elastic 
modulus of the system, in particular in mixtures contain- 



ing a smaller volume fraction of small spheres, x s < 0.5. 
In addition, we compare the dynamics of a sample under 
shear to its quiescent state, showing that the driving intro- 
duced by shear leads to an acceleration of the non-affine 
particle motions, inducing glass melting. A discussion of 
the non-linear rheology of these mixtures and compari- 
son to predictions of mode-coupling theory are reported 
in separate publications [11, 12]. 

METHODS 
Samples 

Suspensions of poly-methylmethacrylate (PMMA) 
particles sterically stabilized with a layer of polyhydrox- 
ystearic acid (PHS) were prepared in a solvent mixture 
of cycloheptyl bromide (CHB) and cis-decalin, closely 
matching the density and refractive index of the colloids. 
In the CHB/decalin solvent mixture, the spheres acquire 
a small charge which is screened by adding 4 mM 
tetrabutylammoniumchloride [13]. This system shows 
almost hard-sphere behaviour, with the volume fraction 
= (47C/3 )nR? being the only thermodynamic control 
parameter, with n the number density of particles and 
R the sphere radius. Binary colloidal mixtures with 
8 ~ 0.2, fixed total volume fraction <j) ~ 0.61 and dif- 
ferent mixing ratios were prepared starting from one 
component stock suspensions. The stock suspensions 
were obtained by diluting a sediment of large particles 
of mean size di = 1.76 ±0.02 /xm (relative polydis- 
persity a = 0.057), or small particles of mean size 
d s = 0.350 ±0.004 jum (a = 0.150). The large particles 
were fluorescently labeled with nitrobenzoxadiazole 
(NBD). For the two one-component colloidal stock 
suspensions, the values of the radius and polydispersity 
were determined from the angular dependence of the 
scattered intensity and diffusion coefficient obtained by 
means of static and dynamic light scattering, respec- 
tively, on a very dilute colloidal suspension (0 ~ 10~ 4 ). 
The volume fraction of the sediment of large spheres was 
experimentally determined as follows: A first guess for 
the volume fraction (Prcp of the sediment was obtained 
using simulation results [14]. The sediment was then 
diluted to a nominal volume fraction ~ 0.4 and ob- 
served using confocal microscopy. The imaged volume 
was partitioned into Voronoi cells and the mean size of 
the Voronoi volume per particle calculated. The ratio of 
the particle volume to the mean Voronoi volume serves 
as an estimate of the volume fraction of the sample. 
This was found to be = 0.43 which corresponds to 
§ 1 rcp ~ 0-68. The small spheres were too small to be 
imaged. Therefore their volume fraction was adjusted 
to give an equivalent rheological response to the large 



spheres. For ideal hard spheres, the energy density scales 
as nkgT, so that the shear moduli must be equal in 
these units. The volume fraction of the one-component 
small particles suspension was adjusted to obtain the 
same normalised shear moduli as for the one-component 
large particles. Although their linear viscoelasticities 
are thus within experimental resolution the same, their 
volume fractions could be slightly different, since the 
samples have different polydispersities. Accordingly, for 
intermediate x s , shear moduli are reported in reduced 
units of the energy density. 



Confocal microscopy 

Quiescent State 

Confocal microscopy experiments on quiescent sam- 
ples were performed using a VT-Eye confocal unit 
(Visitech International), mounted on a Nikon Ti-U in- 
verted microscope with a lOOx Nikon Plan-Apo VC oil- 
immersion objective, and a laser with X = 488 nm. Sam- 
ples were contained in vials where the bottom was cut 
and replaced by a coverslip to allow for imaging [15]. 
Stacks of images of 512x512 pixels, corresponding to 
an x-y plane size of approx. 50 x 50 jiim 2 were acquired. 
Each stack was composed of 101 images obtained ev- 
ery 0.2 fim in z-direction, leading to an imaged volume 
of approximately 50 x 50 x 20 jlim 3 per stack. The time 
needed to acquire one stack was approximately 3.8 s. 
Stacks were acquired at a depth of approx. 30jUm from 
the coverslip. Typically for each sample 7 different vol- 
umes were imaged for 1200s during which 300 stacks 
were collected for each volume to follow the dynamics 
of the samples. The stacks were analysed using standard 
routines [16] to extract particle coordinates and trajec- 
tories. Figure 1 shows typical two-dimensional images 
corresponding to a plane in a stack, acquired for samples 
with different mixing ratios x s . 

Under Shear 

Under shear, samples were imaged using a custom- 
built rotational shear cell (a modified version of the 
model described in [17]), mounted on a Zeiss Axiovert 
M200 microscope with a 63x Zeiss Plan Neo Fluar 
water-glycerol immersion objective and equipped with 
a VT-Infinity confocal unit (Visitech International). A 
glass coverslip serves as bottom plate of the cell to al- 
low for imaging with high numerical aperture objectives. 
The glass surface was covered with polydisperse col- 
loidal hard spheres with a size comparable to the large 
spheres to minimise the effects of wall slip [18]. The top 




FIGURE 1. Typical confocal microscopy images of quies- 
cent samples showing the large fluorescently-labelled particles 
only. The total volume fraction of the samples is (j) =0.61 and 
the relative volume fraction of small particles x s are as indi- 
cated. 



of the cell is formed by a metal cone with 14 mm diam- 
eter and 2° cone angle. Images are acquired at a radial 
distance of 7 mm from the center. The plate and cone ro- 
tate in opposite directions, giving rise to a zero-velocity 
plane in the sample, the depth of which can be adjusted 
through the relative speed of the cone and plate. Images 
were acquired with an Andor iXon 897 EMCCD camera 
for 300s, at an average rate of 10 frames per second. Sol- 
vent evaporation was minimised using a solvent sealing 
at the top of the cell. 



Rheology 

Rheology measurements were performed with a 
AR2000ex stress-controlled rheometer, using a cone- 
plate geometry with 20 mm diameter, 2° cone angle and 
0.054 mm gap. A solvent trap was used to minimise 
solvent evaporation during the measurements. The tem- 
perature was set to 20 °C and controlled within ±0.1 
°C via a standard Peltier plate. The effects of sample 
loading and aging were minimized by performing a 
standard rejuvenation procedure before each test: di- 
rectly after loading, we performed a dynamic strain 



FIGURE 2. Radial distribution function g(r) of the large 
spheres, determined using confocal microscopy. The total vol- 
ume fraction of the samples is = 0.61 and the relative volume 
fractions of small particles are x$ = 0, 0. 1, 0.3, 0.5, 0.7, 0.9 (top 
to bottom). Data are shifted along the vertical axis for clarity. 



sweep, i.e. applied oscillatory shear to the samples 
with a frequency co = 1 rad/s and an increasing strain 
amplitude until the sample was flowing. Before each 
measurement, flow of the sample was induced applying 
oscillatory shear at strain y = 300%. Shear was applied 
for the time needed to achieve a steady-state response, 
i.e. the storage modulus G' and the loss modulus G" 
become time-independent, typically 200 s. Successively, 
the linear viscoelastic moduli were measured at 0.1% 
< y < 0.5% (depending on sample) as a function of time 
to monitor reformation of structure, until the moduli 
reached a time-independent value, typically after 100 s to 
900 s (depending on sample). After this, the experiment 
was started immediately. 

RESULTS AND DISCUSSION 
Quiescent Structure 

To understand mixing-induced changes on the cage 
structure of a one-component glass, we used confocal 
microscopy to determine the radial distribution functions 
g(r) of large spheres, the only species which is fluores- 
cently labeled and therefore visible (Figs. 1 and 2). 

The g(r) for x s = is typical of a glass-forming one- 
component suspension with size polydispersity. It shows 
a pronounced correlation peak at r » di, correspond- 
ing to the highest probability of finding particles in the 
first-neighbour shell, and additional peaks at larger r re- 
lated to particles in the successive neighbour shells. For 
a small volume fraction of small particles (x s = 0.1, 
(j>i = (1 — x s )<j> = 0.549) these features remain, but in 
addition a small shoulder to the right of the first max- 
imum is observed. This indicates a perturbation of the 



cage formed by the large spheres. When increasing x s to 
0.3 (0/ = 0.427), the height of the first-neighbour peak 
decreases, which indicates dilution (also evident in Fig. 
1), and that some particles formerly constituting the cage 
are located at larger distances. These particles are found 
at distances d\ + d s (where the shoulder was observed 
at x s = 0.1) and c/, +2d s , as seen from the correspond- 
ing peaks in g(r). This implies that small particles are 
located in between large particles and hence loosen the 
cage structure. In line with this observation the layering 
of large spheres only extends to the third neighbour shell. 
At x s = 0.5 ((j>i = 0.305) particles are mostly located at 
distance c// + d s and also the probability of finding parti- 
cles at di + 2d s is increased. Moreover, additional peaks 
at d[ +nd s are visible. This indicates that at x s = 0.5 a 
first neighbour shell of large spheres does no longer sur- 
round large particles (Fig. 1), and a transition to a cage 
of small spheres takes place. This is consistent with the 
following geometrical argument: Each small sphere of 
radius R s projects on a sphere of radius R/ = R s /S an 
angle 9 = 2arcsin(l /(l + 1/ 5)). The maximum packing 
of small spheres having the centers separated by this an- 
gular distance, i.e. covering the surface area of a large 
sphere, can be calculated as N — 120 [19]. At x s = 0.5 
the number fraction of small particles for each big parti- 
cle is ^s/^i = x s /8 3 (l —x s ) = 125, i.e. on average each 
large particle is covered by small particles for x s = 0.5 
and hence the first neighbour shell and cage of large par- 
ticles disappear. At x s — 0.7 (0, = 0.183) correlations at 
distances c// + n d s dominate and layering beyond the sec- 
ond neighbour shell vanishes due to the pronounced di- 
lution of the large spheres (Fig. 1). Correlations are fur- 
ther reduced atx s = 0.9 ((pi = 0.061) due to the increased 
dilution of the large spheres (Fig. 1). From the evolu- 
tion of the radial distribution function with increasing x s 
one can therefore conclude that the small spheres occupy 
an increasingly larger fraction of the free volume in be- 
tween the large spheres, inducing a distortion of the cage 
of large spheres until a transition to a system dominated 
by the cage of small spheres is observed. 



Quiescent Dynamics 

To explore the effect of the structural changes dis- 
cussed in the previous section on the microscopic dy- 
namics, we investigate the mean squared displacements 
(MSDs) of the large particles, drf, as a function of 
x s (Figure 3). The system of only large spheres (x s = 
0.0) presents an MSD which, within the accessible time 
range, shows no long-time diffusion, i.e. glassy dynam- 
ics. Moreover, the plateau of the MSD corresponds to 
a localisation of the particles on distances of the order 
Vi(x$ — 0) s=s 0.1c/;, which is typical for a cage in a one- 
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FIGURE 3. Mean squared displacements 8rf of the large 
spheres, normalised by their diameter squared (df), determined 
using confocal microscopy. The total volume fraction of the 
samples is <j) — 0.61 and the relative volume fractions of small 
particles are x s = 0.0 (o), 0. 1 (□), 0.3 (A), 0.5 (V), 0.7 (O), 0.9 
(<). 



component glass. The time-dependence of the dynam- 
ics is similar for x s =0.1, but the localisation length is 
slightly larger. This reflects the small perturbation of the 
cage structure (Figure 2). For.*^ = 0.3 a significant accel- 
eration of the dynamics is observed, for times t > 10 s the 
particles are no longer localised and the MSD increases 
sub-linearly with t. It is expected that diffusive dynam- 
ics is established beyond the accessible time scale. The 
acceleration of the dynamics is related to the consider- 
able distortion of the cage of large spheres induced by 
the presence of the small spheres, which increase the mo- 
bility. A comparable time-dependence of the MSD is ob- 
tained at x s = 0.5, but the displacements are smaller, in- 
dicating a stronger localisation. A stronger localisation 
can be associated with the transition to the cage structure 
of small spheres, as also evidenced in the g(r). Note that 
the caging of the small spheres is apparently incomplete 
yet, and therefore the large particles are not localised. 
For the two largest values of x s the large spheres are 
localised by a cage of small spheres. Accordingly their 
MSD again show no diffusion. The localisation length 
is of the order of V[(x s > 0.7,0.9) 0.02c/,, i.e. about 
8vi(xs — 0), which indicates that the particles are indeed 
localised on the length scale of the cage of small parti- 
cles. Note that the plateau values of the MSDs of these 
samples approach the resolution limit of the setup. 



Linear Viscoelastic Moduli 

In order to establish a link between the microscopic 
structure and dynamics of the samples and their me- 
chanical response, we measured the frequency depen- 
dent linear viscoelastic moduli of the mixtures (Figure 
4). The moduli are reported in units of energy density 
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FIGURE 4. Storage modulus G' (full symbols) and loss 
modulus G" (open symbols), in units proportional to the energy 
density {kgT / (R 1 )), as a function of oscillatory Peclet number 
Pea, for samples with <j) = 0.61 and relative volume fractions 
of small particles x s = 0.0 (o), 0.1 (□), 0.3 (O), 0.5 (V), 0.7 
(O), 0.9 (<), 1.0 (A). 

(nksT) ~ kgT I (R 3 ), with n particle density and: 

This representation removes the trivial effect of differ- 
ent average particle sizes for different values of x s on 
the absolute values of the shear moduli. The data are 
shown as a function of the oscillatory Peclet number 
Pe m = TB/T ffl = (6nrj(o(R i ))/kBT which represents the 
ratio between the period of oscillation, T m = \/(0, and 
the Brownian time, x B = (R 2 /Dq), where Do is the free 
diffusion coefficient. 

At large values of Pe ffl , for all samples G" is larger than 
G' . This response can be associated to the in-cage dy- 
namics, i.e. the short time diffusion of a particle in its 
cage. In contrast at smaller frequencies, i.e. longer times, 
the structural relaxation associated with long-time diffu- 
sion allows us to distinguish the response of a glass from 
that of a fluid. 

The one-component systems (x s = 0.0, 1.0) show the re- 
sponse of a glass. The storage modulus G' is larger than 
the loss modulus G" and no crossing of the two mod- 
uli can be observed at low Pe m , indicating that no struc- 
tural relaxation is observed in the accessible frequency 
window. At x s =0.1 the Pe m (frequency) dependence 
of the viscoelastic moduli is similar to that of the one- 
component systems, i.e. still characteristic of a glass, but 
the viscoelastic moduli are reduced by more than an or- 
der of magnitude, despite the only limited structural de- 
formation of the cage induced by the presence of the 
small spheres. This is however consistent with the larger 
localisation length observed in the dynamics, which in- 
dicates a looser cage structure. At x s = 0.3 the reduc- 
tion of the shear moduli is even more pronounced than at 
x s =0.1. Moreover G' and G" become similar, indicating 
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FIGURE 5. (a) Height of the first, g lmax (black, full sym- 
bols), and second, g2max = g(4 + ^s) ( re d, open symbols), 
peaks of g(r) (Figure 2), (b) Localisation length V/ in units 
of d;, estimated from MSDs at t = 4.5s (Figure 3) (c) Storage 
modulus G' (Pe m = 0. 1 ) (black, full symbols, left axis) and ratio 
G'/G" at the same Pe m (red, open symbols, right axis), in units 
of energy density kgT / (P?)) and (d) theoretical prediction [23] 
for changes in <Prcp at 5 = 0. 175 as a function of x s . 



a weaker solid-like response. This is consistent with the 
large structural distortion of the cage manifested in the 
radial distribution function (Figure 2) and with the faster 
dynamics (Figure 3). Compared to x s = 0.3, at x s = 0.5 
the moduli at large Pe m are bigger, but smaller at low 
Pe®, which is due to the stronger frequency dependence 
of the shear moduli. In addition, G' and G" have almost 
identical values. Such a response is similar to that ob- 
served for depletion driven colloidal gels in the vicinity 
of the gelation boundary [20, 21, 22]. Further increas- 
ing x s to 0.7 and 0.9, the moduli become larger and for 
x s = 0.9 approach the values of the one-component sys- 
tems. This is consistent with a transition to a system 
dominated by cages formed by the small particles and 
with structure and dynamics of these samples. Note that 
the residual distortion induced by the large spheres leads 
to a reduction of the overall elastic response of these 
samples. 

The trends discussed above are summarized by plotting 
G' and the ratio G' / G" as a function of x s , at a fixed value 
of Pe m = 0.1 (Figure 5c). The ratio G'/G" attains the 



smallest value at x s = 0.5, which could reflect a transi- 
tion from a system dominated by cages formed by large 
spheres to a system dominated by cages formed by small 
spheres. This interpretation is supported by the trends of 
the heights of the first and second peaks of g(r) (Figure 
5a): Between x s = 0.3 andjc s = 0.5 the first peak strongly 
drops and then remains nearly constant for larger x s , in- 
dicating the disappearance of the first neighbor shell of 
large particles, i.e the large spheres cage. The second 
peak reaches its maximum at x s = 0.5, corresponding to 
formation on average of a shell of small particles around 
each large particle, and then decreases for larger values 
of x s , due to the further intercalation of small spheres in 
between two large spheres, which leads to caging of the 
small spheres. In contrast, the elasticity of the samples, 
represented by G', reaches a minimum at x s = 0.3. This 
could be explained by the larger localisation length of 
the large spheres at x„ — 0.3 (Figure 5b). Furthermore, 
changes in G 1 are considerably larger in systems with 
a larger volume fraction of large spheres. This can be 
rationalized by considering the effects of the inclusion 
of the second component on the structure of the sys- 
tem in the two cases: In systems at small x s , the small 
spheres can be packed in the free volume in between the 
large spheres, including the free volume within the cages. 
This deforms the cage and shifts random close packing 
(Figure 5d, data for 8 = 0.175 taken from [23]). On the 
other hand, addition of large spheres to a system of small 
spheres only affects the order beyond the first shell, i.e. 
beyond the cage, since the large spheres cannot fill the 
space in between the small spheres. This results in struc- 
tural heterogeneity rather than cage deformation, and in a 
small shift of random close packing (Figure 5d, x s > 0.5). 



Dynamics under shear 

We investigated the effect of shear on the motions of 
large particles in a sample with a major relative volume 
fraction of small spheres (x s — 0.9). The mean squared 
displacements 8rj of large particles were determined in 
the quiescent and steady state of shear, for two different 
shear rates y. For the applied shear rates, the time scale 
introduced by shear, 1 /y, is considerably longer than the 
Brownian time Tb associated with the short-time diffu- 
sion of both large and small spheres. This is quantified 
through the Peclet number Pe-j, = (67Ct]y{R 3 ))/k B T. Both 
time scales are smaller than the structural relaxation time 
of the system, which diverges, although activated pro- 
cesses typically lead to diffusion at long times [24]. 
The velocity profiles obtained in the steady state of shear 
are shown in figure 6. They were obtained by determin- 
ing the velocity of the particles from their trajectories. 
The zero-velocity plane is located at about 15 jJ.m into 
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FIGURE 6. Average velocity of the large spheres as a func- 
tion of position z in the gap, for sample with a relative volume 
fraction of small spheres x s = 0.9 and Pet = 4.2 x 10~' (main 

plot) and 6.7 x 10~ 2 (inset). Arrows indicate the location of the 
plane where the dynamics were measured. 



the sample for both shear rates. The velocity profiles 
show a larger velocity gradient below the zero-velocity 
plane than above. Within each band though the velocity 
profile is linear, indicating laminar flow. For the higher 
shear rate the slower band corresponds to y ~ 0.06 s -1 , 
while the faster band to y ~ 0.25 s -1 . The weighted av- 
erage y~ 0.117 s _1 agrees within uncertainties with the 
expected value of y ~ 0. 120 s _1 . Similarly for the slower 
shear rate the slower band corresponds to y ~ 0.0095 s~ 1 
and the faster one to y ~ 0.041 s , with an average of 
0.017 s~ 1 . The formation of the two bands might be due 
to the different roughness of the particles coated surface 
of the bottom plate and the smooth metallic surface of 
the cone. The dynamics under shear were determined 
in a velocity-vorticity plane within the slower bands for 
both shear rates, corresponding for the large spheres to 
Pe'^4.2 x 10" 1 and 6.7 x 10~ 2 , and at about z = 27 /im 
in the sample (arrows in Figure 6). 
The results of measurements with the shear cell setup 
are reported in Figure 7 as MSDs vs. strain J—jt. The 
strain axis for the quiescent MSD was obtained using 
the faster shear rate y ~ 0.06 s . The quiescent dynam- 
ics show a time dependence similar to that obtained us- 
ing the other confocal microscope setup (compare fig- 
ures 7 and 3): particles are localised on the experimen- 
tally accessible time window and no long-time diffusion 
is observed. One can observe though that the localisa- 
tion length is larger for the measurements with the shear 
cell setup. This might be attributed to the combination of 
two factors [25]: the larger noise level of the multi-beam 
VT-infinity confocal microscope, which arises from the 
cross-talk of the fluorescence emission from many dif- 
ferent particles simultaneously excited; the smaller mag- 
nification (63x instead of lOOx) and the larger pixel size 
(0.25 pLva compared to 0.115 jJ.m), which increase the 
uncertainty in the determination of particle coordinates. 



Note also that, in order to compare to measurements un- 
der shear, the quiescent MSDs are measured in a two- 
dimensional plane instead of a three dimensional volume 
as in the other setup. Application of a slow shear rate, 
corresponding to Pe^ = 6.7 x 10~ 2 , induces a significant 
acceleration of the non-affine dynamics of the large par- 
ticles, as shown in figure 7: The particles are initially lo- 
calised on the same length scale as in the quiescent state 
but become delocalised at y > 6 %, with the MSD in- 
creasing first sub-linearly and then linearly with y over 
the remaining range of measured times. The final linear 
increase of 5r 2 sinf {y°^ t) indicates diffusive behavior. 
At the larger shear rate (Pe'y = 4.2 x 10 _1 ) the particle 
dynamics first show localisation on a length scale smaller 
than in the quiescent state and for y > 3 % the MSD in- 
creases more than linearly with time t and might at larger 
y tend to normal diffusion. For the smaller shear rate, the 
cage-deformation introduced by shear enables the ini- 
tially caged particles to diffuse, resulting in the observed 
acceleration of the average single-particle dynamics and 
glass melting. The larger shear rate is sufficiently large 
to possibly induce cage constriction at short times result- 
ing in the lower localisation length of the MSD. More- 
over, the observed super-diffusive behavior could result 
from the transition from highly constrained in-cage mo- 
tions to out-of cage shear induced diffusive motions. The 
observed behavior is similar to the one which occurs in 
one-component colloidal glasses and dense fluids under 
application of a constant shear rate, as shown in experi- 
ments [26, 27, 28, 29, 30], simulations [31, 27, 30] and 
Mode-Coupling theory [27, 30]. In particular, a link be- 
tween shear-induced cage break up and acceleration of 
the dynamics has been found [29]. Upon application of 
shear, the cage increasingly deforms, until the maximum 
elastically sustainable deformation is achieved, where a 
stress overshoot is observed in rheology, and the cage 
opens, allowing for diffusion in the steady state of shear 
where residual structural anisotropy is observed. Before 
steady state is achieved, super-diffusion is observed at 
the transition from caging to diffusion, corresponding to 
cage yielding. When the shear rate becomes sufficiently 
large, cage constriction is continuously induced by shear 
and a super-diffusive regime is observed in the steady 
state [29], similar to what is observed in the mixture 
for the larger shear rate. Note that for x s = 0.9 the cage 
being deformed is that composed of small spheres. The 
Pe S y = S^Pe'y values for the small spheres are 5.4 x 10~ 4 

and 3.4 xl0~ 3 for the slower and faster shear rates re- 
spectively. Cage constriction effects are typically ob- 
served for Pej > 0.1 in one-component glasses [29]. 
For glasses composed by only one species of particles 
the long-time diffusion coefficient Dl in the steady state 
of shear is dominated by the time scale introduced by the 
shear rate and scales as y 8 [26, 31]. The scaling clearly 
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FIGURE 7. Mean squared displacements 8rf of the large 
spheres as a function of strain y = ft, in units of the squared 
large spheres diameter df, determined by confocal microscopy 
for sample with a relative volume fraction of small spheres 
x s = 0.9, in the quiescent state (□) and in the steady state of 
shear at Pe'^ = 6.7 x 1(T 2 (o) and 4.2 x 10" 1 (A). The MSDs 
under shear only contain non-affine particle motions. The strain 
axis for the quiescent state was calculated using the lower shear 
rate y~ 0.06 s _1 . 



does not hold in this case, since the values on the x-axis 
scale with y. The ratio between the Dl values for the two 
applied shear rates apparently scales with a larger expo- 
nent of approximately 1.6, which could be related to the 
peculiar properties of the mixture. 

CONCLUSIONS 

We presented experimental results on the structure, dy- 
namics and viscoelasticity of glasses formed by binary 
colloidal mixtures with size ratio 8 = 0.2 and different 
mixing ratios. Changes in the properties of these glasses 
as a function of mixing ratio can be rationalized in terms 
of a transition from caging of the large spheres to caging 
of the small spheres. In comparison to a glass com- 
posed of only large spheres, mixing a large fraction of 
large spheres with a small fraction of small spheres in- 
duces pronounced changes in the glass state. The cage of 
large spheres is deformed due to the inclusion of small 
spheres in the free volume between the large particles. 
This loosening of the cage results in increased mobil- 
ity of the large particles and an acceleration of their dy- 
namics. Correspondingly a strong decrease of the elastic 
modulus is observed. Further increasing the fraction of 
small spheres, the cage distortion increases as more and 
more small particles fill the free volume. This is consis- 
tent with random close packing occurring at a larger total 
volume fraction [23]. It also results in a further speed- 
ing up of the dynamics and reduction of the elastic mod- 
ulus. At x s — 0.5 on average each large sphere can be 
completely covered by small spheres and leads to a dis- 



ruption of the cage structure of the large spheres. Con- 
comitantly the response of the system starts to be dom- 
inated by caging of the small spheres. This is seen as a 
tighter localisation of the large spheres and a modulus 
which starts to increase again. This trend continues with 
increasing x s . In systems dominated by the cage of the 
small particles, the large spheres reduce the order on the 
intermediate length scale beyond the first shell, i.e. the 
cage. If shear is imposed on a mixture where caging by 
the small component dominates the response, the initially 
frozen dynamics become diffusive in the experimental 
time-window at small shear rates, and super-diffusive at 
larger shear rates. A stronger localisation at short times 
is also observed at larger shear rates. This indicates that 
application of shear induces melting of the glass by fa- 
cilitating out-of-cage diffusion through elongation and 
deformation of the cage, and cage constriction at large 
shear rates, similar to recent results on one-component 
glasses [29, 30]. 
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